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Abstract 

In this paper we consider the classification of dispersive linearizable partial difference 
equations defined on a quad-graph by the multiple scale reduction around their harmonic 
solution. We show that the A\, A-i and A3 linearizability conditions restrain the number of 
the parameters which enter into the equation. A subclass of the equations which pass the 
A3 C— integrability conditions can be linearized by a Mobius transformation. 



1 Introduction 

Calogero in 1991 [2j introduced the notion of S and C integrable equations to characterize 
those nonlinear partial differential equations which are solvable through an Inverse Scattering 
Transform or linearizable through a change of variables. Using the multiscale reductive technique 
he was able to show that the Non Linear Schrodinger Equation (NLSE) 

idtu = K2[u\ = d xx u + p2\u\ 2 u, u = u(x,t), (1) 

appears as a universal equation governing the evolution of slowly varying packets of quasi- 
monochromatic waves in weakly nonlinear media featuring dispersion. The necessary conditions 
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for the 5-integrability is that P2 is real. If, however, the equation is linearizable then p2 must 
be null, the equation has to be linear or linearizable, as the Eckhaus equation [^ll20j . 

By going to a higher order in the expansion, the multiscale techniques have been used to 
find new S— integrable Partial Differential Equations (PDE's) and to prove the integrability of 
new nonlinear equations [7,8, 13J. Probably the most important example of such nonlinear PDE 
is the Degasperis-Procesi equation [6]. The application of this test to the case of linearizable 
equations has not been done yet. However Calogero and collaborators constructed many in- 
teresting linearizable nonlinear evolution equations by applying complicated transformations to 
linear equations [HIE] - 

In the case of discrete equations it has been shown [Tjl9l llll[T6l419p 24] that a similar situation 
is also true. One can present the equivalent of the Calogero-Eckhaus theorem stating that a non- 
linear dispersive partial difference equation will not be S— integrable if its multiscale expansion 
on analytic functions will not give rise to an integrable NLSE at the lowest order. Moreover 
it was shown on examples that a nonlinear partial difference equation will be C— integrable 
if its multiscale expansion on analytic functions will give rise to a linear partial differential 
equation [23] . 

As was shown in [10 H 1H . IT7]. the introduction of multiple scales on a lattice reduces the 
given discrete equation either to a local nonlinear partial difference equation, by imposing a 
slow-varying condition, or to a partial differential equation of infinite order when dealing with 
analytic functions. Here we choose the second alternative because, as was shown in [10|lll|. only 
in this case the integrability conditions of the discrete equation are preserved. 

In this paper we provide necessary conditions for the linearizability of a class of real difference 
equations in the variable u : 1? — > R and its three nearest-neighbors defined on a 1? square- 
lattice 



Q(u n ,m, Un+l,m, Un,m+1 , Mn+l,m+l ) Pi, P%i •••) — 0, 



(2) 



where /3j's are real parameters. Linearizability conditions will be determined through the multi- 
ple scale perturbative development. 



* , n,m+l 



^"n+l,m 



Figure 1: The Z 2 square-lattice where the equation Q = is defined 



We will suppose that eq. ([2]) is linear-affine in every variable, and our classification will be 
carried out up to the simultaneous Mobius transformations 



i-> (Au nt7n + B)f (Cu n>m + D) 



(3) 
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In Section 2 we will briefly discuss the multiple scale expansion for eq. ([2]) and present the 
integrability conditions which ensure that the given equation is a C— integrable equation of j 
order, with j = 1,2,3, i.e. it is such that asymptotically satisfies the C— integrability conditions 
up to third order. In Section 3 we apply them to the eq. ([2]) and present a sequence of theorems 
which give the conditions on the constants under which the system is Ai, A2 or ^3 integrable. 
At the end we present some conclusive remarks. 



2 Expansion of real dispersive partial difference equations 

We now briefly illustrate all the ingredients of the reductive perturbative technique necessary to 
treat difference equations, as presented in [9lll0j. We will consider here the multiscale expansion 
in the case of analytic functions so as to preserve the linearizability of the given equation. 

Introducing multiple lattice scales, under some obvious hypotheses on the C(°°) property of 
the function u n , m and on the radius of convergence of its Taylor expansion for all the n and m 
shifts involved in the difference equation ([2]), we can write a series representation of the shifted 
values of u n , m around the point (n, m). Choosing 

rii = e ni n, e n% = A^e 4 , 1 < i < K n , rrij = £ mj m, e mj = Mje j , 1 < j < K m , 

where the various constants Aj, Mj and e are all real numbers and we will assume K n = 1 and 
K m = K (eventually K = +00), the total shift operators T n , T m can be rewritten in terms of 
partial shift operators T n , Tn as 

( e ) + °° N j 

3=0 J ' 

K +00 

I'm = Tm J^J Trrij J = Tm ^ ^ ^ i (^^) 
3=1 3=0 

K +00 

T t -TT T (£ni ] TT T 3 - T T V ^ (a c ) 

± n ± m — I n I m >n\ II I mj — In • m / ^ c ^n,mi V*^ 1 ) 
3=1 3=0 

and in terms of differential operators A m \ An) m . These operators are given by appropriate 
combinations of ^jrd mk (explicit expressions and formulae can be found in 0, see [22] for more 
details). We can assume for the function u n m 

= u(n,m,ni,{mj} . 1 ,e) a double expansion in 

harmonics and in the perturbative parameter e 



+00 7 

u n , m = E E £lu i ] ( n i> m 3>3 > l)e w ^ hn -^ m \ (5) 

7=1 0=-7 



with u\ (m,mj,j > 1) = u\ (ni,mj,j > 1) in order to ensure the reality of u n ,m- Then, 
inserting the explicit expressions of the shift operators in terms of the derivatives with 

respect to the slow variables, eq. ([2]) turns out to be a PDE of infinite order. Moreover we 

(0) (0) 
require the functions u\ to satisfy the asymptotic conditions lim ni _ >= |- 00 u\ = 0, V7 and 9, and 

the index 7 chosen > 1 in order to let any nonlinear part of eq. ([2|) to enter as a perturbation 

in the multiscale expansion. 
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2.1 The orders beyond the Schodinger equation and the C-integrability con- 
ditions 

The multiple scale expansion of a nonlinear partial difference equation on analytic functions 
will thus give rise to continuous PDE's. So a multiple scale linearizability test will require that 
an equation is C— integrable if its multiple scale expansion will go into the hierarchy of the 
Schrodinger equation. To show this we need to consider the orders beyond that at which one 
obtains for the harmonic the Schrodinger equation. References [7lll4^fT5] contain different 
approaches to the higher order multiscale expansion of S— integrable nonlinear PDE's. 

Following [7J we can remove all secular terms from the reduced equations order by order and 
thus, in agreement with Degasperis and Procesi [8], we can state the following theorem proved 
by Michela Procesi |21| : 

Theorem 1. A nonlinear dispersive partial difference equation is C— integrable only if its mul- 
tiscale expansion is given by a uniform asymptotic series such that 

1. The amplitude evolves at the slow-times m a , a > 2 according to the a—th equation of 
the linear Schrodinger hierarchy; 

d m JV = (-i)^B a 8° ni uM, B. = -l^, (6) 
where B a are constants. 

2. The amplitudes of the higher perturbations of the first harmonic Uj, j ' > 2 evolve at the 
slow-times m a , a > 2 according to certain linear, nonhomogeneous equations when taking 
into account the proper asymptotic boundary conditions. 

dmjp - {-it-^B a d° ni uf = M a uf> = (7) 

V j, a > 2, where Ba-d^u^ is the a—th flow in the linear Schrodinger hierarchy All 
the other v*p , > 2 are expressed in terms of differential monomials of u'p, p < j. 

If the asymptotic expansion is, besides a necessary condition for integrability, also a sufficient 
condition, has been only conjectured in [21j . The proof of the theorem takes into account that 
for a C— integrable equation by definition we have a linearizing transformation. The next step 
is to assume that the function satisfying the linearized equation admits a similar perturbative 
expansion as in ([5]) and then to perform a multiscale analysis of the linearizing transformation. 

In eq. (J7J) f a (j) is a nonhomogeneous nonlinear forcing term and B a , a > 1, introduced in 
eq. (J6]) are complex constants. Eq. Q represents a hierarchy of compatible evolutions for the 
function . It is obvious that the operators M a defined in eq. ([7]) commute among themselves. 
Thus, once we fix the index j > 2 in the set of eqs. (|7|), their compatibilities imply the following 
conditions 

M a f a ,(j) = M a ,f a (j), Vo-,a'>2, (8) 

where, as f a (j) and f a i (j) are functions of the different perturbations of the fundamental har- 
monic up to degree j — 1, the time derivatives dt a , dt a , of those harmonics appearing respectively 
in M a and M a i have to be eliminated using the evolution equations (JU [7J up to the index j — I. 
The commutativity conditions ([8|) turn out to be an integrability test. 

To construct the functions /<f we define homogeneous spaces of functions of Uj, its n\- 
derivative and its complex conjugate. 
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Definition 2.1. A differential monomial p it*- , j > 1 in the functions Uj\ its complex 
conjugate and its n\- derivatives is a monomial of "gauge" 1 if it possesses the transformation 



property 

P 



4" 



Definition 2.2. ^4 finite dimensional vector space V v , v > 2 is the set of all differential poly- 
nomials in the functions it,- , j > 1, their complex conjugates and their n±- derivatives of order 
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v in e and gauge 1 such that 



order ( d%, u^) = order (d% iA ) = p + j, p>0; 



Definition 2.3. V u {p), p > 1 v >2 is the subspace ofV v whose elements are differential 
polynomials in the functions u^s, their complex conjugates and their n\- derivatives of order v 
in e and gauge 1 for 1 < j < p. 

The basis monomials of the spaces V u (p) can be found, for example, in |22j . 

The functions f a {j) belong to a well defined finite dimensional space T > u {p) and if the relations 
([8]) are satisfied up to the index k, k > 2, we say that our equation is asymptotically C— integrable 
of degree k or A^ C— integrable. 

2.2 C— integrability conditions for the Schrodinger hierarchy 

In this subsection we present the conditions which must be satisfied for the asymptotic C— integrability 
of order k or the A^ C— integrability conditions with k = 1,2,3. To simplify the notation, we 
will use for vtp the concise form u(j). 

The A\ C— integrability condition is given by the absence of the coefficient p2 of the 
nonlinear term in the NLSE. 

The A2 C— integrability condition is obtained by choosing j = 2 in the compatibility 
conditions ([8]) with a = 2 and a' = 3: 

M 2 / 3 (j) = M 3 f 2 (j) , (9) 

where M 2 = d m2 + iB 2 dl 1 and M 3 = d m:j + B^d^. 

In this case / 2 (2) 6 Va(1) and / 3 (2) G V 5 (l) with dim(P 4 (l)) = 2 and dirndl)) = 5, so 
that / 2 (2) and / 3 (2) will be respectively identified by 2 and 5 complex constants 

/ a (2) = au ni (l)|n(l)| 2 + 6u ni (l)u(l) 2 , (10a) 
/ 3 (2) = a|n(l)| 4 U (l) +/3K 1 (l)| 2 u(l) + 7 n ni (l) 2 u(l) + (10b) 
(l)u(l) 2 + e|u(l)| 2 u nini (l). 

In this way, eliminating from eq. ([9]) the derivatives of n(l) with respect to the slow-times m 2 
and m 3 using the evolutions ([6]) with a = 2, 3 and equating term by term, we obtain that 
the A 2 C— integrability conditions give no constraint for the two coefficients a and b. As a 
consequence one can say that the A\ C— integrability condition /9 2 = automatically implies A2 
C— integrability, or that it indeed represents an A2 C— integrability condition. The expression 
of a, /3, 7, 5, e in terms of a and b are: 

a = 0, /3 = ^— , 7 = — , (5 = 0, e = 7. (11) 
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where, for convenience, from now on we will write B for S3. 

The ^3 C— integrability conditions are derived in a similar way setting j = 3 in eq. ©. In 
this case we have that / 2 (3) G V 5 {2) and / 3 (3) G V G {2) with dim(7> 5 ( 2 )) = 12 and dim(:P 6 (2)) = 
26, so that /2(3) and /s(3) will be respectively identified by 12 and 26 complex constants 

/ 2 (3) = n |u(l)| 4 «(l) + r 2 |n ni (l)| 2 n(l) + r 3 |n(l)| 2 (1) + r 4 u (i) u (iy 

+ r 5 u ni (l) 2 u(l) + r 6 u ni (2)\u(l)\ 2 + T 7 u ni (2)u(l) 2 + r 8 u(2) 2 u(l) + r 9 |n(2)| 2 (12a) 
+ ri u(2)n ni (l)u(l) + T llU (2)u ni (l)u(l) + ri 2 n(2)u ni (l)u(l), 

/3(3)=7i|n(l)|V(l)+72|«(l)| 2 n(l) 2 ^n 1 (l)+73k(l)| 2 (1) 

(1) +75|^n 1 (l)| 2 n„ 1 (l) +76^ 

+ liu nini (l)u ni (l)u(l) + 7 8 tin im (1)^ + 79|-"(1)| 4 ^(2) 

+ 7l0 k(l)| 2 n(l) 2 n(2) + 7ll u ni (l)u(2) 2 + 112 u ni (l)\u(2)\ 2 + 7l3 \u ni (l)| 2 n(2) 
+ 7i 4 |u(2)| 2 n(2) + 7l5"Uni (l) 2 u(2)+ 7l6 |«(l)| 2 (2) + 7i 7 t t (l) 2 (2) 
+ 7i 8 n(2)u 

nini (l)«(l)+7i9«(2)ti nini (l)u(l)+ 720«(2)« nini (1)«(1) 
+ 7 2in(2)u ni (2)«(1) + 722^(2)n ni (2)u(l) + 7 23 n nii (2)u ni (l)fi(l) 
+ 724«m (2)n ni (l)-u(l) + j2BU ni i(2)u ni (l)u(l) + 7 26 u ni (2)ii(2)u(l). 

Let us eliminate from eq. ([9]) with j = 3 the derivatives of u(l) with respect to the slow-times 
U12 and m3 using the evolutions ([6]) respectively with o = 2, 3 and the same derivatives of u(2) 
using the evolutions ([7]) with a = 2, 3. Equating the remaining terms term by term, the ^4 3 
C— integrability conditions turn out to be: 

n = -775- [b (rn - 2r 6 ) + ar 7 ] , 6r 7 = - (6 - a) (m + no - r 6 ) + ar 7 , 
4i>2 2 

(13) 

ar 8 = 6r 8 = 0, ar 9 = 6r 9 = 0, ar 12 = a (no - r u ) + 6t 6 + ar 7 , v 7 

(6 - a) T12 = (6 - a) n . 

Sometimes a and 6 turn out to be both real. In this case the conditions given in eqs. (|13p become: 

(b - a) (Rn + R w - R 6 - 2R 7 ) = 0, (6 - a) (I n + I w - J 6 - 2/ 7 ) = 0, 
(6 - a) #8 = 0, (6 - a) J 8 = 0, (6 - a) Rg = 0, (6 - a) Ig = 0, 
a (i? 12 + R u - R 10 - R 7 ) = bR 6 , a (J 12 + hi - ho - I 7 ) = bl 6 , 
{b - a) (R 12 - Rio) = 0, {b - a) (J 12 - J 10 ) = 0, 
where tj = Rj + ilj for j = 1, . . . , 12. The expressions of the 77 as functions of the Tj are: 



(14) 
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7i = (are - 4S.B 2 ti + br 12 ) , 72 = ^§ (K + 0^7) , 

3iB 3 r 3 3iB 3 r 2 
73 = -^^, 74 = 0, 75 = -^, 

3i-B 3 r 4 3i5 3 r 5 
76 = 5 — , 77 = 75 , 78 = 73 5 — , 79 = 0, 710 = 0, 

-t>2 £>2 

3i5 3 r 9 3iB 3 Tn 3i-B 3 Ti 2 

711 = 0, 712 = 2Br' 713 = 2Jh~ ' 714 = 715 = 25a - ' 

3i-B 3 r 6 3i-B 3 Ti 3iB 3 r 8 

7i6 = ' 717 = 718 = ' 719 = 2B2 — ' 720 = 7l5 ' 721 = I^T" ' 

3iB 3 Tj 

722 = 712 , 723 = 716 + 719 , 724 = 713 , 725 = „ ' 726 = 0. 

-D2 

The conditions given in eqs. (|13l appear to be new. Their importance resides in the fact 
that a C— integrable equation must satisfy those conditions. 



3 Dispersive affine-linear equations on the square lattice 

In this Section we derive the necessary conditions for the linearizability of a dispersive and 
homogeneous affine-linear equation defined on the square lattice and belonging to the class ([2]). 
The most general multilinear equation of the class ([2]) has at most quartic nonlinearity. Let us 
introduce into its linear part the solution u n ^ m = K n Q m where K = e 1K and £1 = e~ lul ( K \ We 
get that this equation is dispersive if homogeneous and written as 

Q± = fl iK,m ± «n+l,m+l) + °2 

+ (ax - a 2 ) U^mUn+i^m + (ct\ + a 2 ) U n>m+ iU n+ i im+ i 
+ (Pi - P2) Un,mU n ,m+l + (Pi + P2) 'Un+l,m'Un+l,m+l 
+ llU n ,mUn+l,m+l + 72^n+l,m^n ,771+1 

(15) 

"I" (£l ?3) "n,m%+l,m"n,m+l ~\~ (ft ~\~ ft) W-n,m'Un+l,m^' , n+l,m+l 
+ (ft - £4) «ii+l,mMn,m+l w n+l,m+l + (£2 + £4) ^n,m^n,m+lWn+l,m+l 
C^n,m^n+l,m^n,m+l^n+l,m+l — 0, 

where a\,a 2 £ M\ {0}, |ai| 7^ 1 0,2 1 , are the coefficients appearing in the linear part while 
ai, a 2 , fa, fa, 71, 72, £1, £4, C are eleven real parameters to be determined by using the multi- 
scale procedure described in Section 2. The linear dispersion relation is given by 



u(k) = arctan 



(a\ — a 2 ) sin k 
a i + a i) cos K + 2aia2 



(16) 



Let us look for those transformations which leave the class of equations Q± invariant. They 
will give the equivalence conditions for our classification. It is well known that polynomial equa- 
tions are invariant under simultaneous Mobius transformations u n ^ m \-t (Au n ^ m + B)j (Cu n ^ m + 
D). However, our class of equations is homogeneous with a restriction on the coefficients of the 
linear part. After a Mobius transformation a constant term 

a = B A C + 2B 3 D (ft + ft) + B 2 D 2 [ 7l + 72 + 2 (a x + ft)] + 2BD 3 (a x + a 2 ) (17) 
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0L\ + Q-2 



U n +l,m+l 




Figure 2: Representation of the quadratic nonlinearities of Q± 



will appear, and thus, if we do not want to restrict the coefficients of the equation, we have to 
set B = to have oo = 0. Under the Mobius transformation with B = the coefficients of Q± 
become 



a\ i->- D 3 ai, a 2 ^-D 3 a 2 , ct\ >->■ D 2 [a\ + C (a\ + a 2 )] , a 2 >->■ -D 2 a 2 , (18) 

fa ^ D 2 [ft + C (oi + a 2 )] , ft -> £> 2 ft, (19) 

7l h+ D 2 ( 7l + 2Cai) , 72 ^ D 2 ( 72 + 2Ca 2 ) , (20) 

6 ' y DCi + \CD [3C (ai + a 2 ) + 7i + 72 + 2 (ai - a 2 + ft)] , (21) 

6 ^ £>6 + \CD [3C (ai + a 2 ) + 71 + 72 + 2 («i + a 2 + ft)] , (22) 

6 £>6 + ^C-D [C(ai - 02) + 71 - 72 + 2ft] , (23) 

£4 -> ££4 + [C(ai - 02) + 71 - 72 - 2ft] , (24) 

C ^ C + C 2 [2C (ax + o 2 ) + 71 + 72 + 2 («i + ft)] + 2C + £2) . (25) 



So our equivalence transformation with respect to which we will be classifying the equation Q± 
is a restricted Mobius transformation of the form u n ^ m H> it n>m / '{Cu n%m + D). 

In this article we limit ourselves to consider the case of Q + . The case Q_ will be dealt with 
in a subsequent publication. 

Imposing that p 2 = we get he following Proposition: 

Proposition 1. The lowest order necessary conditions for the linearizability of equations Q+ 
give 6 different classes of equations characterized by different non-superimposed ranges of values 
of the coefficients of the equation Q + . They are: 

• Case LI: 

OL2 = ft = 0, OL\ = fa, 71 + 72 = 2fa, 
t _ t _ (ai+a2) 2 7i72 + (3ai-2a2)a27lft- a i( 2a l-3 a 2)72,3l 
< Si — ?2 — 4aia 2 (ai+a 2 ) ' (26) 

_ _ -(ai-a2)(ai+a2) 2 7i72-a2(-a?-5a 2 ai+2a|)7i/3i+ai(2a2-5a2ai-a|)72/3i 
^3 ?4 4aia2(ai+a2) 2 
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Case L2: 



fh = 0, a\ = Pi, (3ai - 2a 2 ) 0^71 + a\ (2ai - 3a 2 ) 72 = 4ai (a\ - a 2 ) a^Pi, 

(ai+a2)(a|7j-af7|)+2a 2 (af-2a|)7i/3i+2ai(2aj-a|)72/3i-6ai(ai-a2)a 2 /3i 

4aia2(a 2 — a?,) 
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2aia2(ai+a2)(7i -72)^1 + (ai-a2)(a27i- a l72) 2 +2aia2(a 2 -ai)/3i 



Case L3: 



4aia 2 (ai+a2) 




0, ai = 2a 2 , 

602 



71 = 
£3 



e 



(ai+/3i), 72 = g (ai + /Si) 
(ai-2/3i)(2ai-/3i) 
18a 2 



(27) 



(28) 



• Case Z^: 



"2 



6 = 



2a\ = a 2 , 

-5ai/3 1 +a 2 +/3j ! 
6ai 



7i = 
£3 



i(ai+/?i), 72 = |(ai+A) 

(ai-2/3i)(2ai-/3i) 



1801 



(29) 



• Case L5: 



• Case 



"2 

6 
. £3 

a 2 = 
£1 = 
£3 = 



/3 2 , ai = 

3^-2fe^i+y3| 
6ai ' 

18ai ' 



-/3 2 , ai = /3i. 

3/3g+2/3 2 /3i+/3| 
6a 2 

18a 2 



2ai = a 2 

£4 



71 



Mi 

3„ ' 



72 



3/3 2 +2/3 2 /3i+/3 2 
6ai ' 
/3i 2 +6/3 2 /3i+7/3| 
18ai 



Mi 

3 ' 



«i = 2a 2 , 71 = 

? 4 ~~ 18a 2 



2/3i 
72 = ^, 



(30) 



(31) 



The corresponding six subclasses of equations are invariant under restricted Mobius transforma- 
tions of the form u n 1— > u n /(Cu n + D). 

Proof: Following the procedure described in Section 2 we expand the fields appearing in equa- 
tion Q_|_ using the definitions (0 U]) . The lowest order necessary conditions for the integrability 
of Q+ are obtained at the order e 3 of the multiple scale expansion. At this order we get the 
m 2 -evolution equation for the harmonic ujj 1 ' , that is in general a nonlinear Shrodinger equation 
of the form 



id m2 



(i) 



Bod, 



n 2 U l 



(1) 



P2U 



(i)uM|2 



0, where, 



n 2 = n\ 



dui 
dn 



mi, 



(32) 



where the coefficients B2 and p 2 will depend on the parameters of the equation Q + and on the 
wave parameters k and oj, with u expressed in terms of n through the dispersion relation (I16p , 
Let us give just a sketch of the construction of eq. (I32p . omitting all intermediate formulas. 
At 0{e) we get: 



for 9 = 1 a linear equation which is identically satisfied by the dispersion relation (|16p . 
for 6 = a linear equation whose solution is Uj° = 0. 
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At 0(e 2 ), taking into account the dispersion relation (|16|) . we get: 

• for 9 = 2 an algebraic relation between and 

• for 9 = 1 a linear wave equation for , whose solution is given by Ui [n%, m\ ,7712) = 
U\\"n2, fn^), with ni given by eq. ([32]) . 

• for = an algebraic relation between and u\ ■ 

Let us stress here that at C(e 2 ) we find that all the harmonics will depend on the slow- variables 
n\ and m\ through n-i- 

At C(e 3 ), for 9 = 1, by using the previous results, one gets the NLSE (|32p with 

aiq 2 (aj - opsin w 
(af + + zaia2 cos kJ z 



where 



sin k 



7^(°) + 7^W cos K + Tif) cos 2 K + ^j 3 ) cos 3 K + n (A) CQS 4 K 



(01 + a 2 )(a 2 + a| + 2a\a,2 cos k) 2 [(ai — a2) 2 + 2aia2 cos + cos k)] ' 



(33) 



^ _ TZf ] + cos k + n { 2 ] cos 2 « + TZ { 2 ] cos 3 K + TZ^ cos 4 « + ^^ 5) cos 5 K 
(01 + a2)(a 2 + a| + 2aia 2 cos k) 2 [(ai — CI2) 2 + 2aia2 cos k(1 + cos n)] 

(i) (i) 

Here the coefficients 1Z\ , < i < 4, and 7£ 2 , < z < 5, are polynomials depending on the 
coefficients ai, a 2 , ai, a 2 , /3i, /3 2 > 71)72, £i,---,£4 an d their expressions are cumbersome, so that 
we omit them. 

Note that B2 is a real coefficient depending only on the parameters of the linear part of Q+, 
while P2 is a complex one. Hence the linearizability of the NLS equation (132h is equivalent to 
the request P2 = V k, that is 

7e| ) = 0, < i < 5, j = 1,2. (35) 

Eq. (I35p is a nonlinear algebraic system of eleven equations in twelve unknowns, the coefficient £ 
not appearing at this order of the multiple scale expansion. By solving it with the help of the 
computer algebra software Mathematica one gets the solutions (f26H3~l"j) . These solutions are 
computed taking into account that 01,02 E R \ {0} with |ai| 7^ 1 0,2 1 • Let us point out that the 
equations have been solved using only rational algebra, avoiding the unreliability of computer 
algebra calculations when using irrational functions. 

(i) 

One first solves the six equations corresponding to 7£ 2 = 0, < i < 5. Two of them can be 
solved for £1 and £3 in (rational) terms of the remaining ten coefficients. The resulting system 
of equations turns out to be ^-independent and linear in the four variables a±, /3i, 71 and 72. 
Therefore we may write the remaining four equations as a matrix equation with coefficients 
nonlinearly depending on ot2, P2, a\ and 0,2- The rank of the matrix is three. Solutions (I26H31P 
are obtained by requiring that the matrix be of rank 3, 2, 1 and 0. In this way we get 

• Case 1: 

at2 = fa = 0, , , 

4l — 42, S3 — ?4- 
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Case 2: 



OL2 = @2, &1 = Pi, 

ai = 2a 2 , 
7i = 272, 
, - 62) = -ai(6 



(37) 



-2a 2 72- 



Case 3: 



Case 4: 



«2 = -/?2, «i = Pi, 
a 2 = 2ai, 
72 = 271, 

> 01(6 - 6) = - £4) = -«27i- 



a\ = P\ = |(1 + ai/a 2 )72, 
a27i = «i72, 

- £2) = -a 2 7i, 
, ai(& - £4) = /?27i- 



(38) 



(39) 



Case 5: 



Case 6: 



(02 - oi)/3 2 = (02 + ai)"2, 
2aia 2 (ai - a 2 )ai = (ai + a 2 )(72af - 
2a!0 2 /3i = 71 «2 + 72a?, 
(o 2 - ai)(£i - £ 2 ) = (71 - 72)02, 
(a 2 - ai) 2 (^ 3 - £4) = [7 2 (a 2 - 3ai) 

(a 2 + ai)/3 2 = (a 2 - ai)a 2 , 

2aia 2 ai = 71 a 2 + 72a?, 

2aia 2 (ai - a 2 )/?i = («i + 02X720? 



7i 02) 



(40) 



71(01 - 3o 2 )] a 2 . 



7i al) 



(41) 



(a^ - af)(£i - £ 2 ) = [71 ( a i - 3a 2) - 72(02 - 3oi)] a 2 , 
(01 + o 2 )(6 - £4) = (72 - 7i)«2- 

(1) 

Then we impose the remaining five equations 1Z\ = 0, < i < 4 to each of the six obtained 
solutions ([36H4T]) . By a straightforward computer aided computation we conclude that Case 
1 has four subcases, Cases (f26H29|) . that pass those conditions while Case 4 has two subcases, 
(|30l I3T1) . All the other cases do not satisfy the equations = 0, < i < 4. Cases 1 and 4 
represent the necessary and sufficient conditions for the coefficients a and b in (jlOap to be real. 
A direct calculation proves invariance of the resulting equations with respect to the reduced 
Mobius transformation. Let us stress again that Cases (|26H3ip are A2 C— integrable. 

□ 

As a consequence of this proposition we can state the following obvious but important Corollary: 

Corollary 1. If the coefficients a%, o 2 , a\, a%, Pi, P%, 71,72, ,£4 °f equations Q + do not 
satisfy the conditions in Eqs. h2(M31\) then Q + is not linearizable. 

Note that the trivial linearizability condition cx\ = a 2 = Pi = Pi = 7i = 72 = £1 = £2 = £3 = 
£4 = is contained in eqs. (f26H3T|) (cases LI to L6). 

A particularly interesting case is when we consider an equation with at most quadratic non- 
linearity. In this case we get the following Proposition: 
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Proposition 2. // £1 = £2 = £3 = £4 = *n equations Q + then the lowest linearizability 
conditions are: 

• Cases QLla and QLlb: 02 = fa = 0, a\ = fa, 2of — 7a2ai + 2a^ = 44> ai = 
i (7 ± ^33) a 2 , 71 = (l T ft, 72 = (l ± ^/J) ft- 

• Cases QL2a and QL2b: a 2 = fa = 0, a\ = fa, 4af - 12a 2 af + 9a2af - 12a|ai + 4a| = 44> 
ai = (l + 73 ± 2\/* + 3v ^) a2 > 71 = (19 + 18y/2)fa, 72 = ±2^/4 (19 + 18^2) fa. 

Cases QLla, QLlb are in case LI and cases QL2a, QL2b are in case L2. 

Let us consider now the higher order terms of the expansion for the cases L1-L6. Imposing 
the A3 C— integrability conditions (|14p on the real and imaginary parts of the coefficients Tj, 
j = 1, ... ,12, defined in the expression (|12ap . one obtains that the most general A 3 C— integrable 
equation is represented by: 

a (gl + Q 2) 71 a n a 2 7l 
Oil = fa = 7, > «2 = fa = 0, 72 = , 

2«i ai 

6 = 6 = 8^? ' 8^ • 

The case (|42p is the intersection of the cases (|26p and (|27p . As a consequence of the result (|42p 
we have that an equation of the form Q + which satisfies the Proposition i.e. dispersive with 
at most quadratic nonlinearity, will never be C— integrable. 
Moreover it is easy to prove the following theorem: 

Theorem 2. The equations Q+ satisfying the A3 C— integrability ( QSU can be linearized by a 
real Mobius transformation 



7^n,m + & 



if and only if C = - 1 — ^ — • The coefficients of the linearizing transformation are 
Aa'[ 

a n a ^ 

P = ^ 7 = "^7' 

while the resulting linearized equation is 

V n ,m + V n+ i m+ i H {V n +l,m + V n m+l) = 0. (43) 

«1 

Equation (|4"3l) is the most general linear dispersive equation defined on the square. As the 
Mobius transformation used is restricted Mobius transformation, (|43p can be assumed to be the 
canonical equation for this case. 

4 Conclusions 

In this work we presented the complete classification of the linerizable dispersive partial differ- 
ence equations belonging to the Q + class using multiple scales expansions around a periodic 
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discrete wave solution of the linearized equation, up to the fifth order in the perturbation ex- 
pansion parameter. The resulting linearizable system depends on three free parameters, and 
requesting the resulting equation to explicitly linearize via a Mobius transformation reduces 
the free parameters to one. It is still an open problem the proof of the linearizability of the 
equation Q + satisfying just the A 3 C— integrability conditions ([32]) . Maybe by going higher in 
the perturbation expansion we could be able to fix the parameter £ according to the theorem of 
last Section. 

This calculation shows that the multiple scale expansion can be effectively used to classify 
discrete equations. 

Work is in progress for the derivation of the integrable class of dispersive partial difference 
equations belonging to the Q + class. They are provided by a reduced set of equations with 
respect to the one considered in this case, as p2 7^ is just a real constant. Moreover the 
integrability conditions are given by a larger number of nonlinear algebraic equations for the 
parameters entering into the equation and, as such, are much harder to solve. 

Work is also in progress in the study of the Q_ case. In this case we get at the lowest order 
in the perturbative parameter a nonlinear system of partial differential equations relating the 
zeroth and the fundamental harmonic. The solution of this equation is the key ingredient in the 
classification of this class of equations which contains all dispersive equations belonging to the 
ABS classification of multilinear equations on the square. 
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